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O ■ Abstract 

o , 

, We consider a mechanism for the formation of superhumps in the TV Col system, 

P5 ' based on the possible existence of a precessional spiral wave in the accretion disk 

^ ^ . of the system. This mechanism can act in binaries with arbitrary component-mass 

ratios, and our precessional spiral wave model can be applied to explain observed 
superhumps of all types. 
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1 Introduction 



^ I Superhumps are modulations of the light curves of binary systems with periods that 

^ I differ from the orbital periods by several percent, and are observed mainly during su- 
O ■ peroutbursts in SU UMa systems. The main observational features of superhumps are 
described in [1]. Various authors have put forward models to explain the superhump phe- 
nomenon (a brief description of various models and their problems are given, for instance, 
O I in [1]). Currently, the most popular model explains these light variations as an effect of 
' precession of the outer regions of the accretion disk. The presence of the Lindblad 3:1 
resonance in the disk results in an instability that leads to precession of the outer part of 
the disk, with a period that is appreciably longer than the orbital period. The beating of 
^ , the orbital and precessional periods gives rise to the periodic variations that are identified 
^ ! with superhumps. This model has several shortcomings, the most important being that 
it implies a limit on the maximum component mass ratio. In order for the Lindblad 3:1 
resonance to be located inside the accretion disk, the ratio of the donor and accretor 
masses q must be lower than ~ 0.33 [2], while there are several observed systems with 
superhumps that have appreciably higher q values, such as TV Col. 

TV Col is a "permanent superhumper" , or one of a class of systems in which su- 
perhumps are always observed. It was suggested in [3] that cataclysmic variables with 
high mass-transfer rates (TV Col is such a system) may persist in a superhump regime 
that is not interrupted by quiescent periods. Among stars with superhumps, TV Col 
has an unusually long orbital period. Port ~ 5.5 h [4]. Since there is a direct correlation 
between the orbital period of a system and the donor mass (see, for instance, [3]), this 
may be as high as ~ O.56M0, while the component-mass ratio is probably in the range 
q ~ 0.6-^0.9 [5]. It is obvious that the Lindblad 3:1 resonance cannot be located inside the 
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accretion disk with such a q value. This casts doubt on the popular superhump model 
indicated above. 

In 2004, a new mechanism for the formation of supcrhumps in SU UMa systems was 
suggested in [6,7]. The basis of this mechanism is the idea that a precessional-type density 
wave can form in the accretion disk. 

The precessional wave in the accretion disks forms as a result of interactions between 
elliptical streamlines. The asphericity of the gravitational field of the binary system 
forces the semimajor axes of the streamlines to precess opposite to the direction in which 
the matter flows. The rate of this precession is proportional to the semimajor axis of 
the corresponding streamline. Since the streamlines in the disk cannot intersect, their 
interaction results in the establishment of a certain equilibrium precession rate of all the 
streamlines, with their apastrons lining up so that they form a spiral pattern. Because 
the matter velocity is minimum at the apastrons, the density of the matter at these points 
grows, and a precessional spiral density wave is formed in the disk. 

After the formation of the precessional density wave in the disk, the rate of accretion 
grows sharply (by up to an order of magnitude). Matter approaches the surface of the 
accretor along the precessional wave, and the region of accretion is localized in azimuth, 
and, hence, forms a radiating spot at the surface of accretor. The increase in the accretion 
rate due to the density wave explains both the development of a superoutburst and the 
amplitude of the superhump. The wave slowly precesses in the stationary (observer's) 
coordinate frame, leading to a shift of the region of enhanced accretion (i.e., of the 
superhump) with each rotation of the system. The beating of the orbital period and the 
precessional period of the wave result in the superhump period, which is slightly larger 
than the orbital period. 

This superoutburst model based on the presence of a precessional spiral wave in the 
accretion disk [7] made it possible for the first time to explain all important observational 
manifestations of superoutbursts and superhumps in SU UMa systems, including their 
period, duration, and energy, the anticorrelation of the brightness and color temperature 
at the maximum of an ordinary superhump, the late- superhump phenomenon, etc. More- 
over, this model does not place strict constraints on the component mass ratio, and so 
can be applied to the superhumps in the TV Col system. 

The main aim of the present paper is to investigate the possible formation of preces- 
sional spiral waves in accretion disks in systems with q > 0.33. We will also consider the 
possibility of explaining all types of superhumps with this model. 

2 The Model 

The model used for our numerical simulations of mass flows in binary systems is 
described in [8]. The flow is described by a three-dimensional system of gravitational 
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gas-dynamics equations, including nonadiabatic radiative heating and cooling of the gas: 



< 




(1) 



-pv grad $ + p2m;2 (r(T, T^,) - A(T)) . 



Here, p is the density, v = (m, t>, w) the velocity vector, P the pressure, e the internal 
energy, $ the Roche potential, nip the proton mass, flbm the orbital angular velocity of the 
system, and r(T, T^^) and A(T) the radiative heating and cooling functions. The system 
of gas-dynamical equations was closed with the ideal gas equation of state, P = (7 — l)pe, 
where 7 is the adiabatic index, which was specified to be Y3. 

This system of equations was solved numerically using the Roe-Osher method [9- 
11] adapted for a multiprocessor computer. A three-dimensional, noninertial, Cartesian 
coordinate system corotating with the binary was used for the modeling. 

Only half of the space occupied by the disk was modeled, due to the symmetry of the 
problem about the equatorial plane. The computational domain was specified to contain 
both the disk and stream of matter from Li, and had a size of 1.12^4 x 1.14^4 x 0.17A 
(where A is the binary separation). The computational grid contained 121 x 121 x 62 
cells and was distributed among 81 processors in a twodimensional (9 x 9) matrix. The 
grid became more dense closer to the equatorial plane, making it possible to resolve the 
vertical structure of the disk. 

The boundary conditions at the outer edges of the computational domain were spec- 
ified as follows. In all cells except those in the vicinity of Li, we assumed a constant 
density, pb = lO^^pLi, where pi^ is the density at Li , a temperature 13600° K, and zero 
velocity. The accretor was specified as a hemisphere with a radius of 10~^A. Matter 
entering the cells forming the accretor was taken to fall onto the star. The stream was 
defined as a boundary condition: matter with a temperature of 5800° K and an x ve- 
locity Vx = 6.3 km/s was injected into a region with radius 0.014A centered on Li. The 
density of the matter at Li was specified so that the mass-transfer rate corresponded to 
the observed rate in TV Col. 

We made three computational runs with different parameters of the system. Table ^ 
gives the accretor (Mi) and donor (M2) masses and the binary separation (A) for the 
three models. The binary separation was adjusted to be consistent with the observed 
orbital period of the TV Col system of 5.5 h. Since the component mass ratio for TV Col 
is not known precisely, we used the two extreme values q = 0.93 and q = 0.3 (which 



No. 



1 
2 
3 



0.3 1.4 0.42 1.923 

0.7 0.8 0.56 1.745 

0.93 0.6 0.56 1.655 



Table 1: Model parameters q, Mi, M2, A 
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0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.2 1.3 1.4 1.5 
X/A 



0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.2 1.3 1.4 
X/A 



Figure 1: Density contours and velocity vectors in the equatorial plane of the system 
in models with q = 0.7 (left) and q = 0.93 (right). The accretor is located in the center 
of the region shown in the plot. In the vicinity of the precessional wave, the density 
contours are shown using a gray scale. The maximum density corresponds to black. 

still allow the disk to contain the 3:1 Lindblad resonance), as well as the intermediate 
value q = 0.7. For the given rate of mass inflow into the system, the accretion rate in 
the model corresponded to ~ 10~^ Mq. The modeling covered large intervals of time (up 
to several dozen orbital periods), making it possible to reach a quasistationary accretion 
regime with a constant mass of the disk. 

3 Results of The Modeling 

Figure H shows density contours and velocity vectors in the equatorial plane for model 
nos. 2 and 3 from Tabled 

Places where the density contours become closer together correspond to stationary 
shocks in the disk. A dense, round accretion disk and a compact circumdisk halo form in 
the system. At the edge of the circumdisk halo, which is formed by matter that collides 
with the stream when the former moves around the accretor, the density rapidly falls to 
the background values. The interaction of the circumdisk halo gas and the stream results 
in the formation of a shock (hot line; see [6-9,12-16]) outside the disk. A twoarmed tidal 
spiral shock forms in the disk, with both the arms located in the outer regions of the disk. 
A precessional spiral density shock is clearly visible in the inner regions of the disk. The 
two-armed spiral shock is at rest in the noninertial coordinate frame fixed to the binary, 
while it moves with a slow velocity in the inertial observer's frame. The period of the 
precessional wave in the noninertial coordinate frame is slightly longer than the orbital 
period. 

To find the rate of precession, we determined the phases at which the maximum den- 
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Figure 2: Dependence of the angular velocity of the density peaks at three different 
distances from the accretor in models with q = 0.7 (left) and q = 0.93 (right). Angular 
velocities are shown for a noninertial coordinate frame rotating with the orbital angular 
velocity of the system. The horizontal dashed line corresponds to the orbital velocity of 
the binary. 



sity is observed for various distances from the accretor. The precession of the wave results 
in a shift of these phases with time. The typical time dependences of the angular displace- 
ments of the density wave for three distances from the accretor (the beginning, middle, 
and end of the wave) are shown in Fig. El These plots demonstrate that the velocities of 
the streamlines in the disk level out, and the wave moves with some equilibrium velocity 
^Iwav in the coordinate frame fixed to the binary. 

Table |21 lists the precession rates ^Iwav found via the numerical modeling for systems 

with various component mass ratios q. Table |21 also presents the period excexxes e"*" = 
p+ p 

^^p — where Pj~^ is the superhump period. 



4 Connection Between Superhumps and Precessional 
Waves 

As a rule, five types of superhumps observed in cataclysmic variables can be distin- 
guished [17]. 

(1) Positive (or ordinary) superhumps have a period that is a few percent longer 
than the orbital period. Such superhumps were first observed in SU UMa during a 



No. 


Q 


^wav / ^bin 


e+ 


1 


0.3 


-0.07 


7% 


2 


0.7 


-0.1 


11% 


3 


0.93 


-0.05 


5% 



Table 2: Rate of precession in the three models 



5 



Figure 3: Formation of a negative superhump. The positions of the precessional wave 
(W) and bright (hot) spot (H) are shown for three different times. The directions in 
which the wave and leading edge of the superhump move are shown by arrows. 

superoutburst. 

(2) Orbital superhumps represent a modulation of luminosity with the orbital period. 

(3) Late superhumps are observed after the end of a superoutburst, and have the same 
velocity as the positive superhump, but shifted by half a period relative to the latter. 

(4) Permanent superhumps have the same features as positive superhumps, but are 
observed in stars that lack supcroutbursts. 

(5) Negative superhumps have a period that is several percent shorter than the or- 
bital period. Such superhumps were first discovered during monitoring of systems with 
permanent superhumps. 

The precessional spiral wave model is able to explain all types of observed superhumps. 
The main observational features of positive and late superhumps are a direct consequence 
of the formation of a precessional spiral wave in the disk, as is explained in [7]. The orbital 
superhump can be explained by the release of energy in stationary shocks in the disk and 
circumdisk halo, such as the hot-line and tidal shocks; the existence of tidal shocks is not 
precluded in the precessional wave model. Permanent superhumps can also be explained, 
if the mass-transfer rate in the system is high enough to sustain a high accretion rate 
due to prolonged existence of the precessional spiral wave. The existence of luminosity 
modulations whose period is shorter than the orbital period, i.e., of negative superhumps, 
can also be plausibly explained in this model. 

In the model, the superhump radiation arises from a relatively compact region near 
the surface of the accretor, into which matter flows along the precessional wave. If the 
bright spot arising in this region is located above the accretor surface (i.e., the rotation 
of accretor does not influence the motion of the spot) and has a tendency to spread 
out due to diffusion, the leading edge of the spot may be observed somewhat earlier 
after each rotation of the system, creating a modulation with a period shorter than the 
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orbital period. The fact that the observed periods of negative superhumps do not display 
significant scatter can be taken to justify the assumption that the energy release occurs 
above the surface of the accretor. 

The necessary rate of diffusional spreading can be estimated using the observed period. 

Let us denote the observed "period excess" for a negative superhump e~ 

where Porb is the orbital period of the system and P~f^ is the observed period of the 
negative superhump. For TV Col, P^^ ~ 5.2 h [4], so e~ ~ —0.055. If the observed 
period is the result that of beating between the orbital period and rotational period of 
the leading edge of the spot Pieadi this period can be found using the formula 

Plead = = —Porb = — ~ 95.3 h. 

-r orb ~ i^sh ^ 

The hnear velocity of the spot edge will be 

27rr* 27rr* e" 
viead = s — = — 5 r ~ 0.18 km/s, 

-Head -Torb 6 + -1- 

if the radius of the accretor is ~ 10^ cm. The velocity of the leading edge has two 
components: the rate of diffusional spreading v^iff and the (negative) velocity of the spot 
center, relative to the retrograde precession of the wave, v^av ■ Knowing the period of 
the precessional wave P^av (or the period excess of the positive superhump e"^), we can 
find the linear velocity of the wave, Vwav'- 

27rr* ^ PX - Porb 27rr* e+ n ^ i / 

v^av = —B— = -27rr.^^- = —— ■ -—- ^ -0.4 km/s. 

^wav ^sh orb ^orb t -p J- 

This velocity is negative, since the direction of the wave precession is retrograde. Ac- 
cordingly, Vdiff will be given by 

Vdiff = Viead - Vwav ~ 0.58 km/s. 

Let find a relation between the diffusion rate derived from observations and the pa- 
rameters of the disk gas. The time scale for diffusion over a distance L is r = L'^/D, 
where D is the diffusion coefficient [18]. If r = Porb, the linear diffusion rate is 



Vdiff = - = 



D 



^ V Porb 

In this case, the diffusion coefficient for the TV Col system will be 

^ 1.1 • 10^^ f e- e+ V 1^ 9 , 

The diffusion coefficient for an ionized gas is defined as D — I • Vrms, and depends on 

on 1 n6 rp2 

the mean free path / and rms velocity Kms- For ionized gas, I — ^ l^A ' where T 

1.3 • lO^T^/^ 



is the temperature of the gas, n the number density of particles, and A 
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T(K) 



Figure 4: Dependence D{T) for various n. The dashed hne corresponds to the D observed 



the gas particles. Accordingly, the diffusion coefficient is a function of the temperature 
and number density of the matter: D = f{n,T). Thus, if we know the value of D from 
observations, it is possible to find n and T. 

Figure m shows the dependence D{T) for various n. The dashed line shows the value 
of D that corresponds to the observed diffusion for TV Col. If we can estimate the typical 
temperature from observations of the superhump, Fig.|4]can be used to find the density of 
matter in the spot. Let us now assume that there is pressure equilibrium in the accretion 
disk, so that nT = noTo, where uq and Tq are the number density and temperature of 
the outer parts of the disk. Assuming that Tq ~ lO'^K, we can find the typical number 
density in the disk, no. Figure |S1 shows a set of solutions for various hq (for Tq = lO'^K). 
For convenience, we also show |e~(T)| for TV Col in this figure. If we suppose the 
temperature of the spot to be T ~ 10^ K and use the observed value of |e~| (dashed line), 
then no ~ lO^^cm"^. The estimated values of T and no are within the range typical 
for accretion disks in close binary systems [8]. Thus, this model with the diffusional 
spreading of the spot can adequately explain the observed negative superhumps. 

5 Conclusions 

Our numerical modeling has shown that precessional density waves can form in the 
accretion disks of systems with large component mass ratios, up to g = 0.93. 

Our derived dependence of the precession period on q is complex. Figure IHl shows 
the theoretical dependence e+(g) obtained in [19] (bold dashed line). The open circles 



in TV Col. 




where k is Boltzmann's constant and m is the mass of 
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Figure 5: Theoretical dependence |e (T)| in TV Col for various values of uq. The dotted 
line corresponds to the |e~| observed in TV Col. 

in this figure show the period excesses found for three different sets of parameters of 
TV Col. We also show the model e~^{q) values for OY Car and IP Peg found in [6,7]. 
Despite the significant scatter, all the e+(g) values are fairly close to the theoretical line, 
although the values for q = 0.49 (IP Peg) and q = 0.93 (TV Col) arc outliers. Note that 
the rotational rate of the precessional density wave depends not only on the component 
mass ratio, but also on the size of the region occupied by the wave. The wave can exist 
only in a gas-dynamically unperturbed region of the disk, whose size is limited by the 
depth of penetration of stationary shocks - the two arms of the tidal shocks and the hot 
line (extended region of interaction of the stream from Li and the circumdisk halo) - 
into the disk. The high masstransfer rate in TV Col may result in an increase in the 
extent of the hot line and, consequently, to precession rates that are lower than expected 
theoretically. It is interesting that, in the case of high component mass ratios, the size of 
the unperturbed region may be significantly limited by tidal waves, which would result 
in an even larger decrease of the precession rate with increasing q, compared to the 
theoretical value. Nevertheless, the dependence is close to the theoretical dependence, 
and has a peak in the range ~ 0.4 0.6. 

In TV Col, superhumps with the period Pgh ~ 6.3 h [4] are observed, so that the 
superhump period excess is ~ 15%. This value is out of the range of the theoretical 
estimates of the precession rate (but docs not exceed the value obtained in [6] for IP Peg). 
The theoretical dependence for e+(g) was obtained in the model that did not take into 
account the gas pressure; i.e., it effectively estimates the precession rate of the semimajor 
axis of the orbit of a free particle. The excesses of the observed and simulated precession 
rates over the theoretical values may indicate that the estimated disk radius used in the 
simulations, raise — [2]? is too low (in fact, the derivations in [2] did not take into 
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Figure 6: Dependence of the superhump period excess on the component mass ratio. 
The bold dashed hne corresponds to the theoretical dependence given in [19]. The thin 
dashed line corresponds to the same dependence for a disk whose radius is 10% larger. 
The open circles show the period excess of the precessional density wave we have found 
for various possible values of q for TV Col. The triangle corresponds to the period excess 
found for OY Car in [7], and the square to the period excess found for IP Peg in [6]. 



account the gas pressure). If the disk radius were 10% larger (r^j^c = j^^)) all the 
derived precession rates would be below the theoretical line (shown by the thin dashed 
line in Fig. IH)). 

In summary, we conclude that the presence of a precessional density wave in the disk 
can lead to superhumps both in SU UMa stars and in binaries whose component mass 
ratios prohibit the location of the 3:1 Lindblad resonance inside the disk. The precessional 
wave model is able to explain all types of observed superhumps. 
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